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■ We investigate a flat anisotropic (5+l)-dimensional cosmological model in the presence of 

the Gauss-Bonnet (GB) contribution in addition to usual Einstein term in the action. We 

compared it with (4+l)-dimensional case and found a substantial difference in corresponding 

cosmological dynamics. This difference is manifested in the probability of the model to 

have smooth transition from GB-term-dominated to Einstein-term-dominated phases — this 

' probability in a reasonable measure on the initial condition space is almost zero for (4+1) 

I , case and about 60% for (5+1) case. We discuss this difference as well as some features of 

5h , 

the dynamics of the considered model. 
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12 . 1 131 ? ]) showing a revival of interest to this theory after about a decade of 



3,Q, 



ow progress (for achievements of 80-th and early 90-th of the last century see, for example, 



17| ) • The main feature of this extension of the General Relativity is that the resulting 



OO 

equations of motion are second order differential equations, in contrast to 4-th order equations in 
other modified gravity theories. Another important feature of Lovelock gravity is finite number of 
terms in the Lagrangian in contrast to string gravity, where there is an infinite row of additional 
terms. This means that in perturbative string gravity we can study effects arising from corrections 
to GR only at the level of small perturbations, otherwise it would be necessary to include all infinite 
set of terms into analysis. In Lovelock gravity it is possible to study regimes with non-Einstein 
terms being equally or more important than the Einstein part of the action. 

The Lagrangian of Lovelock gravity in an iV-dimensional space consists of all terms which are 
topological invariants in lower dimensional spaces. As a result, there are no corrections to Einstein 
theory in (3 + 1) dimensions: the second curvature invariant, Gauss-Bonnet term, appears in four 
dimensions and have a contribution to the equations of motion for N > 4. The third invariant, 
cubic in curvature, starts to contribute from (6 + 1) dimensions etc. In the present paper we study 
a flat anisotropic (5 + l)-dimensional Universe, so the Lovelock Lagrangian has only one additional 
term in the form of the Gauss-Bonnet combination. 
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In the multidimensional anisotropic cosmology the GB-term leads to two qualitatively new 



features. A new type of sin, 
derivatives appears Q. 



paper 
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;ularity with finite scale factors and energy density and diverging time 



21] . In the recent 



20l | . and recollaps of a flat Universe becomes possible 
22| it was shown that in the dynamics of a (4 + l)-dimensional anisotropic Universe these 
two regimes are in some sense typical, and in most cases they prevent the Universe from reaching a 
low-energy regime from a standard Big Bang singularity. Only special initial conditions with three 
equal Hubble parameters allow smooth evolution from Big Bang singularity to a 4-dimensional 
Kasner regime (which is exact solution in Einstein gravity), all other trajectories either describe 
recollapsing Universe or meet the nonstandard singularity. 

In the present paper we provide a similar analysis for a flat (5 + l)-dimensional Universe in 
Lovelock gravity. We show that the former result is valid only for that particular number of 
space-time dimensions, and though in the (5 + l)-dimensional case there are no additional terms in 
Lovelock Lagrangian in comparison with the (4 + l)-dimensional one, the smooth evolution from 
high-energy to low-energy regimes in (5 + l)-dimensional Universe is typical and does not require 
severe fine-tuning of the initial conditions. 

We consider a flat anisotropic Universe with the metrics = diag(—l, a 2 ,b 2 , c 2 , d 2 , e 2 ) and the 
action 



where the Gauss-Bonnet term 



S = j ^g(R + aGB)d 6 x, 



(~itd jyiklm p a r>ik jy _i_ r>2 

Ltd — K riiklm — Kik + K . 



Introducing five Hubble parameters in a standard way, it is possible to write down equations of 
motion in the form of first integral 



2H a H b + 2H a H c + 2H a Hd + 2H a Hf + 2Hf,H c + 2HbHd + 2HbHf + 2H c Hd 

+2H c H f + 2H d H } + 24a [ H a H b H c H d + H a H b H c Hf + H a H b H d H f (1) 

+H a H c H d Hf + H b H c H d Hf] = 0, 
and five dynamical equations. The first equation of motion has the form 
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2(H b + H 2 ) + 2(H C + H 2 C ) + 2(H d + Hj) + 2(H f + Hj) + 2H b H c + 2tf 6 tf d 



+2H b H f + 2H c H d + 2^ c fr/ + 2H d H f + +8a 



3H b H c H d Hf + (Hi + H 2 ) 



x (H c H d + H C H/ + Hdfly) + (H c + H 2 c )(H b H d + H b H f + H d H f ) 



(2) 



+(H d + H 2 d )(H b H c + FfeF/ + H c Hf) + (if, + Hf)(H b H c + tf 6 i? d + F c tf d 



23(]. 



0. 



four other equations can be obtained by cyclic transmutation of indices 

Two obvious limiting cases can be identified easily. First, in the low-curvature regime, when 
the contribution from the Gauss-Bonnet term is small, we have well-known generalized Kasner 
solution. In this solution scale factors have power-law dependence on time with the power indexes 
Pi, and, correspondingly, Hi = pi/t with the Kasner conditions Y,Pi = 1 and Y,Pi = 1- 

On the other hand, when Einstein part of equations of motion are negligible in comparison 
with the Gauss-Bonnet contribution, another form of power-law solution exists 23f] . There are two 



form of this solution, both require ^ PiPjPkPl = 0, and either = 3 or ^ PiPjPk = 0. 

i>j>k>l i i>j>k 

The latter case contains solutions with pi = (a, b, 0, 0, 0) where a and b are arbitrary numbers. In 



addition to discussion in [23] we can show now that this form is the only possible solution of the 
second class up to transmutation of power indices. 

Indeed, expressing p\ from Yl PiPjPk = 0, then substituting it into Yl PiPjPkPl = and 

i>j>k i>j>k>l 

expressing p2 we get 



1 P4P5 + P3P4 + P3P5 ± P3P4P5 V -3pjpj ~ 2p3P'jP5 ~ ^PzPjPi ~ Zpjpl ~ ^pjPiPb ~ 3P 3 Pj 

_2^2 9999 99999 

2 P3PIP5 + P 3 P\ + P%PiPb + P3PIPA + PlPl + PlPl 

(3) 

We want to show that the expression under the radical is non-positive. To show this consider the 
case when this expression crosses zero: 



Det = -(3p 4 p 5 + 2p 3 p 4 p 5 + 2p 3 p 5 p 4 + 3p 3 p 4 + 2p 3 p 4 p 5 + 3p 3 p 5 ) = 



(4) 



It would happen when 



P3 



p 4 p^(p 4 + P5 ± 2y /z 2p| 



P4P5 



3pj + 2p 4 p 5 + 3p{ 



(5) 
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One can easily verify that radical expression in ([5]) is negative for any nonzero p^ and p§. This 
means that the expression under radical in (3) has the same sign for any combination of non-zero 
Pi, an easy check shows that this sign is negative, ruling out all these combinations. If one of pi is 
zero, the condition Yl PiPjPkPi = requires that at least one of remaining indexes is also equal 

i>j>k>l 

to zero, and in this case the condition ^2 PiPjPk = further requires that the third index also 

i>j>k 

vanishes. As a result, the only possible form of the studied class is the combination (a, b, 0,0,0) 



found in 



2M 



This means that this class of solutions is a rather special one, and it is not surprising that we 
do not see it in our numerical simulations. From now on we consider only first class of power-law 
solution as a solution with standard cosmological singularity in a high-energy regime. 

Apart from high-curvature and low-curvature asymptotic power-law regimes the system (1-2) 
describes complicated behavior in the range of Hubble parameters where both Gauss-Bonnet and 
Einstein contributions are important. We have studied this system numerically, starting from initial 
conditions distributed randomly in the interval [0,1] for three Hubble parameters, in the interval 
[-1,1] for one of them and calculating fifth Hubble parameter from the constraint equation (1). The 
integration stopped when one of the following regimes are reached: 

• For expanding Universe: low-curvature Kasner regime, nonstandard singularity or recollaps. 

• For contracting Universe: high-curvature power-law regime ("standard singularity") or non- 
standard singularity. 

Matching two branches, expanding and contracting, of numerical solution allows us to construct a 
full history of the modeled Universe evolving from one of listed states to another. 

The main feature which distinguishes the (5+1) case from the (4+1) one is that now the 
transition from a standard singularity to low-curvature Kasner regime does not require severe 
fine-tuning (we remind a reader that in the former case only trajectories with three scale factors 
equal to each other can show this type of transition) . 

In Fig. 1 we presented a typical behavior of the model with "standard singularity" — > Kasner 
smooth transition. In (a) panel we show Hi{t) evolution curves, in (b) panel — ^iij(i); individual 

_ i 

Kasner exponents pi are shown in (c) and their sum ^2pi(t) in (d). One can easily see that this is 

_ i 

a high-energy power-law i^2,Pi{t) = 3) to low-energy Kasner {^2,Pi{t) = 1) transition indeed. We 

i i 

have found that about 60% of all the trajectories belong to this class. Studying examples of this 
kind of dynamics for various initial conditions we also noticed that final values of power indexes 
are not distributed in homogeneous way on the Kasner sphere, but prefer the situation clearly 
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FIG. 1: Typical behavior of the model with smooth high-to-low-energy power-law regimes transition. In (a) 
we present the behavior of individual Hi (t) , in (b) — the expansion rate ^2 Hi {t) , in (c) — individual Kasner 

i 

exponents pi(t),'m (d) — the sum ^2pi{t), and in (e) — the sum PiPjPk- 

i i>j>k 
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FIG. 2: Typical behavior of the model which experiences a non-standard singularity. In panel (a) we 
presented Hubble functions Hi(t) and in (b) — the denominator of Zfj. 

seen in the Fig. 1(c) with one negative index and four positive and close to each other. Inserting 
pi ?s p2 ~ P3 ~ P4 into the Kasner conditions one can found that Yl piPjPk ~ —0.3. This 

i>j>k 

estimation meets the results of our numerical analysis fairy well, which can be seen from Fig. 1(e) 
where the combination PiPjPk is plotted. The reason for pi to be approximately equal in 

i>j>k 

the low-energy stage remains a mystery and may indicates that some hidden symmetries we are 
ignorant about are involved, so the problem requires more detailed investigation. 

One of "non-smooth" solutions indicated above is the case of nonstandard singularity. This is 
the case when sometime during the evolution Hi diverges while Hi remain finite. Solving eqs.(2) 
with respect to Hi we can express them in the form of fraction, and the nonstandard singularity 
occurs when the denominator crosses zero while numerator is regular and non-zero. In Fig. 2 the 
behavior of our model is shown for this case. In panel (a) we present the Hubble functions Hi(t), 
in panel (b) — the behavior of the denominator of the fractional expression for Hi. At t ~ 1.2 
the denominator crosses zero, numerator and Hubble parameters remain regular and non-zero. As 
this singularity is weak according to the Tipler terminology 2j], our numerical program typically 



goes through it, and the numerical results after the singularity are not physically valid: despite 
Hi(t) are regular, Hi(t) are singular, and, hence, scalar curvature invariants diverge, which makes 
the non-standard singularity physical, but not a coordinate one. So, the evolution presented in 
Fig. 2 (a) is only valid till the point t ~ 1.2 where the non-standard singularity occurs. Depending 
on initial conditions, the nonstandard singularity can be found either in future or past history of 
the particular Universe. These trajectories take up to 15% of all the trajectories. 
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Finally, the last case is the recollaps. In this case, starting from expansion (the initial conditions 
in this case are chosen in a way to ensure ^2 Hi > from the beginning), the Universe begins to 

i 

contract at some point and ends up in a singularity. In Fig. 3 we presented a typical behavior of 
the model which experience recollaps. In the panel (a) we showed individual Hi(t), their sum in 
(b), individual Pi(t) in (c). They diverge at some points due to the definition of pf. p% = —Hf /Hi. 
In the panel (d) we presented Hi(t) and one can verify zeros of Hi(t) correspond to "singularities" 
of Pi{t). Also in panel (e) of the Fig. 3 we presented the denominator of the H, so the reader can 
verify that singularities in pi have nothing to do with the nonstandard singularity described above. 
Last thing to mention, this case occupy remaining 25% of all the trajectories. 

We have considered dynamics of a flat anisotropic vacuum (5 + l)-dimensional Universe in 
Gauss-Bonnet gravity. Unlike the (4+ l)-dimensional case studied earlier j^ . we have found that 
nonstandard singularity and possible recollaps which usually prevented five-dimensional Universe 
from reaching a low-curvature regime do not dominate in six dimensions. About 60% of trajec- 
tories in the latter case show smooth transition from GB dominated epoch to a low-curvature 
six-dimensional Kasner regime. The other possibilities are a "non-standard" singularity (about 
15%) and a recollaps (about 25%). As we have used random set of initial conditions, other realiza- 
tions can show slightly different numbers. However, our main result has a qualitative character and 
could not change significantly — a transition from Big Bang to low-curvature asymptotic is rather 
general in (5 + 1) dimensions in contrast to the (4 + l)-dimensional case. This may indicate that 
Lovelock cosmology in the world with odd number of space dimensions is much less pathological 
than in even-dimensional worlds. 

Some dynamical features of the model still require additional investigation. Another possible 
direction of future work is to include a matter into analysis, which can modify the results for both 
(4 + 1)- and (5 + l)-dimensional cases. 

Acknowledgments 

This work is partially supported by RFBR grant 08-02-00923 and scientific school grant 
4899.2008.2. 



[1] D. Lovelock, J. Math. Phys. 12, 498 (1971). 

[2] S. Nojiri, S. Odintsov and S. Ogushi, Int. J. Mod. Phys. A17, 4809 (2002). 
[3] T.Torii and H.Maeda, Phys. Rev. D71:124002 (2005). 



9 



[4] T.Torii and H.Maeda, Phys. Rev. D72:064007 (2005). 

[5] H.Maeda, Phys. Rev. D73:104004 (2006). 

[6] M.Aiello, R.Ferraro and G.Giribet, Class. Quant. Grav. 22, 2579 (2005). 

[7] J. Grain, A.Barrau and P.Kanti, Phys.Rcv. D72: 104016 (2005). 

[8] M.Nozawa and H.Maeda, Class. Quant. Grav. 23, 1779 (2006). 

[9] R. Cai and N.Ohta Phys. Rev. D74:064001 (2006). 

[10] C.Garraffo and G.Giribet, Mod. Phys. Lett. A23, 180 (2008). 

[11] M.Dehghani and N.Farhangkhah, Phys. Rev. D78:064015 (2008). 

[12] M.Dehghani, N.Bostani and S.Hendi, Phys. Rev. D78:064031 (2008). 

[13] A.Paranjape, S.Sarkar and T.Padmanabhan, Phys.Rev. D74:104015 (2006). 

[14] N.Dcruelle, Nucl. Phys. B327, 253 (1989). 

[15] N.Deruelle and J.Madore " Kaluza-Klcin cosmology with the Lovelock Lagrangian" In: Liege 1986 

Proceedings, Origin and early history of the Universe, p. 277-284. 

[16] N.Dcruelle and L.Farina-Busto, Phys. Rev. D41, 3696 (1990). 

[17] T.Kitaura and J.Whccler, Phys.Rev. D48, 667 (1993). 

[18] S.Kawai and J. Soda, Phys. Rev. D59:063506 (1999). 

[19] H.Yajima, K.-I. Macda and H.Okubo, Phys. Rev. D62:034020 (2000). 

[20] A.Toporensky and S.Tsujikawa, Phys. Rev. D65:123509 (2002). 

[21] S.Alexeyev, A.Toporensky and V.Ustiansky, Phys. Lett. B509, 151 (2001). 

[22] R. Chingangban, M. Sami, P. V. Trctyakov, and A. V. Toporensky, Phys. Lett. B661, 162 (2008). 

[23] A. Toporensky and P. Tretyakov, Grav. Cosmol. 13, 207 (2007). 

[24] F.J.Tiplcr, Phys. Lett. A64, 8 (1977). 



